Dispersion of the dielectric permittivity of dense and cold atomic gases 
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On the basis of general theoretical results developed previously in JETP 112, 246 (2011) we 
analyze the atomic polarization created by weak monochromatic light in an optically thick, dense 
and cold atomic ensemble. We show that the amplitude of the polarization averaged over a uniform 
random atomic distribution decreases exponentially beyond the boundary regions. The phase of 
this polarization increases linearly with increasing penetration into the medium. On these grounds, 
we determine numerically the wavelength of the light in the dense atomic medium, its extinction 
coefhcient, and the complex refractive index and dielectric constant of the medium. The dispersion 
of the permittivity is investigated for different atomic densities. It is shown that for dense clouds, 
the real part of the permittivity is negative in some spectral domains. 
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I. INTRODUCTION 

Improvements in techniques for cooling of 
atomic gases in atomic traps make their use very 
promising for practical applications in various ar- 
eas of fundamental science and technology such 
as metrology, development of frequency standards, 
and quantum information problems [H-Q- The 
largest number of applications envisioned for cold 
and ultracold atomic ensembles have at their foun- 
dation the interaction between the medium and 
electromagnetic radiation. This interaction also 
underlies many methods for diagnostics of the 
states of these ensembles. Recently, dense atomic 
clouds, in which the average interatomic distances 
are comparable with the optical wavelength, have 
attracted much attention. This interest is in part 
generated by such fundamental and practically im- 
portant physical effects as Anderson (strong) local- 
ization of light and lasing in a disordered 
medium |14l - [l8| which can take place in dense en- 
sembles. The aim of recent studies in this field 
was to observe these effects experimentally and to 
describe them theoretically. 

One challenging problem in the area of disor- 
dered atomic systems is that the studied atomic 
ensembles normally consist of a large number of 
atoms in samples that are produced with a low 
duty cycle. The larger number of atoms is required 
in order to obtain sufficient signal to noise to study 
the subtle effects of interest. Such experiments 
require realistic modeling in order to extract the 
essential physics of the observed processes. How- 
ever, it is challenging to treat these problems as 
a multi-atom scattering process, and such studies 
have been limited to several thousand atoms [1- 
[2l| : this should be compared with the characteris- 
tic 10^ atom-sized samples of recent experiments. 
Alternative theoretical approaches, even if approx- 



imate in nature, are then desirable. 

The present paper is devoted to the theoreti- 
cal description of optical properties of dense and 
cold atomic clouds. The problem of a dense atomic 
ensemble belongs obviously to the field of macro- 
scopic electrodynamics. The main approach here 
is based on usage of such averaged characteristics 
as the field strength and atomic polarization. The 
key point in a macroscopic approach is in finding 
the susceptibility or dielectric constant of the dense 
ensemble. The influence of density of the medium 
on its susceptibility can be analyzed on the basis 
of the idea of a local field and, following from it, 
the Lorentz-Lorenz formula [22|- This formula is 
sufficient to solve completely the problem of the de- 
pendence on density only if the difference between 
the polarizability of a free atom and its polariz- 
ability in the medium can be neglected [23|. As 
we will show below, for the considered cold atomic 
ensemble this is not the case. The resonant dipole- 
dipole interatomic interaction causes atomic level 
shifts and broadening and thereby essentially mod- 
ifies the atomic polarizability. An explicit analyt- 
ical expression for the susceptibility, which takes 
into account this modification, was obtained ear- 
lier in [2^. The calculation in [2^ was based on 
the relevant macroscopic statistical description of 
the polarization response of the medium to an ex- 
ternal field. Part of the approximations made in 
[20] are valid only for relatively low density en- 
sembles and thus the corresponding results have a 
restricted range of applicability. Constitutive rela- 
tions connecting atomic polarization and an exter- 
nal field can be obtained consistently only in the 
framework of a microscopic approach based on the 
notion of the discrete structure of matter consist- 
ing of separate atoms. 

A consistent microscopic approach has been al- 
ready applied for analysis of infiuence of inter- 
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atomic interactions on spontaneous deca y o f an 
impurity atom embedded in a dielectric [2J-|27|. 
Quite a number of works were devoted to col- 
lective decay in dense homogeneous multiatomic 
media and to properties of spontaneous emission 
of such media initially excited by a weak exter- 
nal field (see, for example [28L'39'| and references 
therein). In these works the main attention was 
focused on the influence of the density of the en- 
semble on its afterglow, i.e. on secondary radia- 
tion. In the present paper we consider the influence 
of interatomic interaction on the properties of the 
ensemble itself. We study the spatial distribution 
of atomic polarization created by weak monochro- 
matic light in a cold atomic ensemble. We show 
that in the case of a uniform random atomic dis- 
tribution amplitude of polarization averaged over 
space configurations decreases exponentially be- 
yond the boundary regions. Its phase increases 
linearly with distance into the medium. On this 
ground we determine numerically the wavelength 
of the light in the dense atomic clouds, its extinc- 
tion coefficient as well as complex index of refrac- 
tion and dielectric constant of the medium. We 
also analyze the dispersion of the permittivity for 
different atomic densities. Note that similar prob- 
lems arise in classical electrodynamics when study- 
ing light scattering from a medium with random di- 
electric inclusions [40]. Such inclusions strongly in- 
fluence light propagation in such media, giving rise 
manifold internal scattering and essentially modi- 
fying its averaged dielectric constant. 

An important feature of the present work is in 
taking into account the polarization properties of 
light. Nearly all the mentioned above papers on 
multiatomic systems used only a two-level model 
for the atoms. This prevents a correct consider- 
ation of the light polarization, adequate descrip- 
tion of the resonance dipole-dipole interaction at 
small distances and, as a consequence, correct cal- 
culation of shifts and broadenings of atomic levels. 
We also do not use a model of an averaged con- 
tinuous medium in our calculations. That is, the 
influence on the intrinsic spatial disorder of the 
atoms in the ensemble is considered. As a speciflc 
illustration of this, recent approaches to atomic 
physics-based localization studies have considered 
systems of reduced dimensionality. One way to 
achieve this for light localization is to optically cre- 
ate a quasi one dimensional system through mod- 
ification of spatially larger samples. Such opti- 
cal channels, with wavelength-scale transverse di- 
mensions, can be created through quantum op- 
tical techniques based on electromagnetically in- 
duced transparency, for instance. Alternatively, a 
strongly focused far-off-resonance laser can gener- 
ate a type of optical wave guide through the dense 
sample, allowing quasi one dimensional localiza- 
tion for a much weaker, but near-resonance probe 



beam. Theoretical modeling the average proper- 
ties of such generated optical wave guides, essential 
to interpretation of experiments, may be done us- 
ing the effective optical responses of the resulting 
medium, as we discuss in the current paper. 

The remainder of this paper is organized as fol- 
lows. In Section 2 we describe our basic physical 
assumptions and the calculational approach. Sec- 
tion 3 presents results of numerical simulations. 
We conclude with a brief synopsis of the results, 
highlighting the main points of the present report. 

II. BASIC ASSUMPTIONS AND 
APPROACHES 

Consider the temporal dynamics of a system 
consisting of iV-fl motionless atoms. Let N atoms 
form the cloud. These atoms are identical and have 
a ground state J — separated by the frequency 
from an excited J = 1 state. The natural line width 
of this state is 7. One atom is located far from the 
cloud and has the same J = O J = 1 struc- 
ture of levels but a different transition frequency 
uJs and a different decay constant 7^. We will as- 
sume that initially all atoms of the cloud are in 
the ground state and the separated atom, which 
we will refer to as a source atom or simply the 
source, is in a coherent state which is a superpo- 
sition of the ground and a small admixture of the 
excited state. In the course of spontaneous de- 
cay such an atom creates an electromagnetic field 
which is a superposition of vacuum and a small ad- 
mixture of a one photon state. As is known, this 
superposition approximates a weak coherent state 
of the field with good accuracy. Under the infiu- 
ence of the field, the atoms of the cloud are excited 
and in due course emit secondary radiation which 
can be absorbed by other atoms of the cloud. The 
process of manyfold photon exchange determines 
the dipole-dipole interatomic interaction and man- 
ifests itself in such phenomena as spontaneous de- 
cay modification, collective atomic state formation 
and so on. 

The microscopic description of dynamics of the 
considered ensemble is based on the non stationary 
Schrodinger equation for the wave function of 
the joint system consisting of atoms and the field 
generated in the process of the evolution 

ih^ = (2.1) 

The Hamiltonian of the system H can be presented 
as the sum of the Hamiltonians of the free atoms 
and the free field Ho, and the operator V of their 
interaction. In the dipole approximation used here, 
we have 

F = -^dWE(ra), (2.2) 
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E(r) = E(+)(r) +E^-^(r 



(-)( 



(2.3) 



. ^ J e^aUka exp(zkr) + h.c. 



V 



where E^^' are the operators of the positive and 
negative frequency components of the field; a^a is 
the photon annihilation operator in a mode with 
wave vector k and polarization a; V is the quanti- 
zation volume; (1*^°^ is the dipole moment operator 
of the atom a, e^a are polarization unit vectors. 

We will seek the wave function ip as an expansion 
in a set of eigenstates {|^)} of the operator Hq: 



(2.4) 



Here, the subscript I defines the state of all atoms 
and the field. 

The key simplification of the approach employed 
is in restriction of total number of states \l) taken 
into account. We will calculate all radiative cor- 
rection up to the second order of the fine structure 
constant. In this case we can consider only the 
following states (see (4l| ) 



\g,g,...g) (g) \ka) 
\g,g, ...g) ® \vac) 
\9,9, ■■■9,e,g, -g) C 
Iff, •■•ff,e,g...g,e,g, 



(2.5) 
(2.6) 

\vac) (2.7) 
.g) ® |ka) (2.8) 



In the rotating wave approximation it is enough 
to take into account only the states (12.51) and 
(j2.7p . States without excitation both in atomic and 
field subsystem (|2.6p allow us to describe coherent 
states of the source atom. Non resonant states 
with two excited atoms and one photon (|2.8p are 
necessary for a correct description of the dipole- 
dipole interaction at short interatomic distances. 
Note that, in considered case, there are three ex- 
cited states for each atom e = |J, m), which dif- 
fer by the value of angular momentum projection 
m = —1, 0, 1. Therefore, the total number of one- 
fold excited states ^^Tf) is 3(iV + 1). 

Equation should be supplemented by an 

initial condition. According to our previous dis- 
cussion, we will consider the case when initially the 
field is in a vacuum state, all atoms of the cloud 
are in the ground state and the source atom which 
we denote by index s is in superposition of ground 
and one of the excited states |J, m). Designating 
the corresponding amplitudes as bp and 6o, we can 
write 



^{O) = h'^\g')+ho\e, 



(2.9) 



where the index Cso corresponds to the one of the 
three possible states of atom s which is populated 
in the initial moment of time. 



In the framework of the assumptions made, the 
amplitude of state -0^/ = \g') does not change dur- 
ing the evolution of the system hgi{t) = fog, because 
transitions to this state from other states taken 
into account are impossible. The transition from 
(|2.6p to any of the states is also impossible. 

To determine all other amplitudes we have to 
solve the set of equations which follows from (|2.ip . 
In spite of the performed restriction of the number 
of states, this set of equations is infinite because 
of the infinite number of degrees of freedom of the 
field. We can however exclude amplitudes of states 
with one photon and obtain a finite closed system 
of equations for be{t) = &e„(t); a ^ s. For Fourier 
components b^{uj) we have (at greater length see 

m) 



■ Wa)(5ee' - Sgg' (l^)] &e' (l^) = Aes(w). 



(2.10) 

Matrix elements See' (w) for e and e' correspond- 
ing to different atoms describe excitation exchange 
between these atoms. Assuming that in state ipf,i 
and ijje atoms h and a are excited correspondingly, 
in the framework of the pole approximation (see 
[43|), we have 



See' (w) 



E 



(2.11) 



3 -Si- 



Here is the projection of the vector r = — rf, 
on the axis of the chosen coordinate system and 
r = |r| is the separation between atoms a and h. 

If e and e' correspond to excited states of one 
atom then See' ('^) differs from zero only for e — e' 
(m ~ ni'). In this case See(i^) determines the 
Lamb shift and the decay constant of correspond- 
ing excited state. Including Lamb shifts in the 
transition frequency uja we get 



See(w) = -i7Q/2- 



(2.12) 



The term Aes(a;) in the right-hand side of Eq. 
(|2.10p describes excitation of the cloud atoms by 
the radiation of the source. Assuming that the size 
of the atomic ensemble is negligible compare with 
the distance from it to the source, and neglecting 
the secondary excitation of the source atom s by 
reradiation from the cloud, we have 



Aes(l^) 
Ses(t^) 



(w - LUs + i7s/2) 
hrs 

X exp {ikrs + ikre) 



Ses(w), 



(2.13) 



(2.14) 
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Here k =a;n/c. Relation (|2.14p is written in a co- 
ordinate frame originating at some point inside the 
cloud; Te are radii locating the atoms; n is a unit 
vector oriented from the source to the cloud. In 
obtaining the expression for I]es(iuj) we took into 
account the non-applicability of the pole approx- 
imation because of the large separation between 
the cloud and source, we used the rotating wave 
approximation for the same reason and kept only 
one term which decreases most slowly with ■ All 
these factors generate the differences between S 
and the elements of matrix S. 

Knowledge of explicit expressions for Aes(aj) and 
See' ('^) allows us to determine the amplitudes of 
all one-fold excited states (|2.7p . Note that system 
(|2.10p can be reduced to an integral equation by 
using the continuous medium approximation. This 
significantly simplifies the solution of the problem 
for a two- level atom system [35l - [39| . Moreover, in 
this case, even an analytic solution is possible for 
spatially homogeneous spherical clouds. This solu- 
tion neglects, however, the important properties of 
real physical systems, and therefore we will solve 
the linear system (j2.10p numerically. In a numer- 
ical solution we can correctly describe all polar- 
ization effects and take into account the random 
inhomogeneities of the medium. 

Introducing the inverse matrix which, as shown 
in [20], is a resolvent operator of the considered 
multi atomic cloud 

i?ee'(w) = [{UJ - UJa)5ee' -See'(w)]"\ (2.15) 

we can write the solution of the system (|2.10p as 
follows 

6e(c^) = ^ee'(w)Ae'.(w). (2.16) 

For amplitude be{t) we get 

oo 6oexp(-iwi) i?ee' (w)Se's(l^) 
h.{t) = I — . 

J 27r (jj -ujs^ ^7s/2 

— oo 

(2.17) 

This relation give us the possibilities to find the 
distribution of excited states at any instant of time. 
In this work we are interested in the spatial distri- 
bution under quasi static conditions. Such condi- 
tions can be realized if decay times of all collective 
states of the dense atomic ensemble are much less 
than the decay time of the source atom s. 

Let us consider the relation (|2.17p for a time 
interval much less that but larger than the 
mentioned collective relaxation times. Formally, 
relations for the quasi steady state regime can 
be obtained by two limiting processes. First we 
should pass to the limit 7s — >■ and then to 



t 00. Realizing these limiting processes and 
taking into account that lim (to — lUs-\- ilsl'^)~^ = 

^(w — Ws), where <j(x) is a singular function 
and that lim <; (cj — ujs) exp(— icji) = —2'r:i5{uj — 

LUs) exp^—iujst) we get 

beit) = boex.p{-iuJst) ^ i?ee'('^s)5^e's(Ws)- 

e' 

(2.18) 

By using be(t) we can obtain amplitudes of all 
states taken into account in our calculations (see 
[19]) and consequently, the wave function of the 
considered system. Among other things, this al- 
lows calculation of the polarization as the aver- 
aged dipole moment of unit volume of the atomic 
ensemble. For a given projection /i of the polariza- 
tion vector we have 

E ('^^^^y (2.19) 

a£AV 

Here is the operator of the corresponding pro- 
jection of the dipole moment of atom a. The 
summation in (|2.19p is made over all atoms lo- 
cated in a mesoscopic volume AV near the point r. 
Quantum-mechanical averaging is performed over 
the wave function of the system. 

In analyzing the polarization it is convenient to 
select positive and negative frequency parts and 
use a basis of circular polarization (fi = 0, ±1): 

V^{r,t)^rl-\r,t)+rl+\r,t). (2.20) 

Using the known wave function and taking into ac- 
count the short life time of the nonresonant virtual 
states with two excited atoms and consequently its 
small contribution to the polarization, we find 

7'(+)(r,t) = p(+)(r)exp(-zc^si); (2-21) 

rl-'\r) = 1^ E E^C"e.K)Se..(...). 

aeAV Cb 

The additional index m at means that under 
summation we have to include only those states 
e™ of atom a which give contributions to the cor- 
responding projection of the polarization vector. 
In the basis of circular polarization such contribu- 
tion comes only from one Zeeman sublevel with 
m — fi. Due to the optical isotropy of the atomic 
ensemble, the orientation of the atomic polariza- 
tion vector coincides with the orientation of the 
polarization of light exciting the atoms. The latter, 
in turn, depends on the specific Zeeman sublevel 
rris of the source atom which was excited initially. 
In the case when the quantization axis coincides 
with the vector n the configurationally averaged 
atomic polarization has only one nonzero projec- 
tion = nis- Thus to determine the polarization 
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we have to take into consideration only the Zeeman 
state with m = rus- 

In the next section, we will use relation (|2.2ip 
to calculate the spatial distribution of atomic po- 
larization and analyze on this foundation coherent 
light propagation through ensembles of different 
densities. 



III. RESULTS AND DISCUSSION 



A. Atomic polarization 



Expression (|2.21l) allows us to consider atomic 
clouds with different shapes and with different 
atomic spatial distributions. We however will fur- 
ther consider mainly model cylindrical clouds with 
random but (on average) uniform atom distribu- 
tions along the vector n. For definiteness let us 
assume that the state TOs = 1 of the source atom 
is excited initially. In such a case the source cre- 
ates a nearly plane right-hand circularly polarized 
wave in the area of the cloud and the vector of the 
atomic polarization has only one nonzero compo- 
nent, which we will refer to without index. 

Fig. 1 shows the spatial dependence of the ab- 
solute value (Fig. la) and phase (Fig. lb) of the 
complex quantity 'P*^+^(r) for different detunings 
A = — of the source probe radiation from ex- 
act bare atomic resonance. The calculations were 
made for a cloud with length L — 10 and radius 
i? = 20. Hereafter in this paper we use the in- 
verse wavenumber of the resonant probe radiation 
in vacuum k^^ — Aa/27r = cjlixLiJa as a unit of 
length. In these units, the mean density of atoms 
is 71 = 0.2. To avoid the influence of boundary ef- 
fects at the lateral surface of the cylinder as well 
as diffraction effects caused by the sharp bound- 
ary we calculate atomic polarization 'p(+) (r) only 
for an area near the axis of the cylinder where we 
can neglect the dependence of the polarization on 
r. In this area we deal with a quasi one-dimension 
case. The polarization depends only on z. Our 
analysis shows that for the considered parameters 
this take place for the inner portion of the cylinder 
with r < 15. Results shown in Fig. 1 are obtained 
by averaging of the atomic polarization over the 
region with radius r = 10. 

The curves in Fig. 1 were obtained by aver- 
aging over the random locations of atoms inside 
the cloud. The total number of statistical tests 
was about 6 • 10^. Despite such a large number 
the curves which were not smoothed additionally 
keep indications of fluctuations. These fluctuations 
manifest themselves most clearly far from the front 
edge of the cloud. Here the averaged polarization 
is extremely small in comparison with the polariza- 
tion corresponding to any random specific spatial 
configuration of the ensemble. 
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FIG. 1: Spatial distribution of atomic polarization, (a) 
Amplitude of polarization abs{T'^^\z)), semi logarith- 
mic scale; (b) Phase of polarization arg{'P^'^\z)). Cal- 
culations were made for a cylindrical cloud with length 
L = 10 and radius R = 20, atomic density is n = 0.2. 



In spite of the fluctuations the results shown in 
Fig. 1 allow us to make several important conclu- 
sions about the spatial dependence of the polariza- 
tion. First, beyond the boundaries and near the 
ends of the cylinder (z = and z = 10) the phase 
of polarization wave increases linearly. Second, 
beyond these areas we have a single-exponential 
decay of the atomic polarization. And last, in 
the boundary regions with size of about 1.5 ^ 2 
we see peculiarities connected with the fact that 
atoms located here interact mainly with atoms sit- 
uated on one side of them, inside the cloud. This 
causes some modification of the dipole-dipole in- 
teratomic interaction. Besides that, the electro- 
magnetic wave reflects from the base edges of the 
cylindrical clouds. This leads to formation of a 
standing wave of polarization. This effect is most 
evident at the far edge of the cloud {z — 10) for 
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the wave strongly detuned from resonance (curves 
corresponded to A = —7 and A = 2.57). For 
these waves absorption is small and the amplitude 
of the reflected wave slightly decreases inside the 
medium. At the very edge of the cloud we have 
either a node or an antinode of the standing wave 
depending on the optical density of the medium. 
In our case there is a vacuum beyond the cloud. Its 
refractive index is equal to unity so for large posi- 
tive detunings we have a node (curve A = 2.57 in 
Fig. la) at the edge and for negative ones we have 
an antinode (A = —7). 

Data shown in Fig. 1 were obtain for several dif- 
ferent detunings and one given density of the cloud. 
Note however that all peculiarities discussed above 
were observed in our calculation for the full range 
of considered parameters; that is, for all different 
detunings and for all considered densities of atoms. 

Knowledge of the polarization of the atomic en- 
semble allows us to make some conclusions about 
light propagation in it. Three averaged quantities, 
these being the polarization, the field strength, and 
the electric displacement are proportional to each 
other. The coefficients of proportionality for re- 
gions away from the boundaries cannot depend on 
the spatial coordinates because here we deal with 
a quasi uniform medium. The linear increasing 
of the phase of the polarization and the single- 
exponential decay of its amplitude means that in 
the corresponding the spatial dependence of 
the polarization and of the averaged field strength 
£ are as follows 



r{z) = Voexp{i{k' + ik")z); 
£{z) = exp(i(fc' -f 



(3.1) 



Here we have taken into account that only one 
component of each vector is nonzero and that these 
components depend only on z. The real k' and 
imaginary k" parts of the wave number can be de- 
termined from the decay coefficient and wavelength 
of the polarization, i.e. from the angles of inclina- 
tion in the region of the linear dependence of the 
curves shown in Fig. 1. The results of correspond- 
ing calculations are depicted in Fig. 2 

Fig. 2 shows how interatomic interactions mod- 
ify the spectral dependencies of absorption and 
reflection in dense media. In dilute media, both 
absorbtion and refraction indexes increase linearly 
with density according to the relation k = ko + 
na{uj)/2, where n is the density, k and (t{uj) are 
the complex wave number and the complex cross 
section of light scattering from free atoms. In the 
case when atomic motion and atomic collisions can 
be neglected the latter gives a Lorentz profile for 
the absorption coefficient and a corresponding dis- 
persion curve for the refractive index. The influ- 
ence of collective effects causes essential distortion 
of the spectra. The absorption spectrum is non- 
symmetric, there are noticeable shifts of the maxi- 
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FIG. 2: Spectrum of imaginary (a) and real (b) parts 
of the wave number of a plane electromagnetic wave 
in atomic ensembles of different density; fco — tda/c is 
the wave number of the resonant source radiation in 
vacuum. 



mum of absorbtion which, in the considered range 
of densities, are in the blue wing. The amplitude 
of the absorbtion increases slowly with density and 
there is an evident tendency towards saturation. A 
density increasing from 7i = 0.2to7i = 0.5 leads 
only to a 25% increase of maximum absorption, 
which is much less that under density increasing 
from n = 0.1 to n — 0.2. Saturation effects in 
our interpretation connects with level shifts caused 
by strong dipole-dipole interaction for dense me- 
dia. These shifts cause also essential nonhomoge- 
neous broadening of the spectral profiles, this being 
clearly seen in Fig. 2. 

Calculation of the real and imaginary parts of 
the wave number permits us to assess qualitatively 
the loffe-Regel criterium for strong light localiza- 
tion in atomic media. According this criterium, lo- 
calization can be observed if the transport length 
of a photon is less than its wave length in the 
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medium. The wave length of the photon is deter- 
mined by the inverse real part of the wave number 
A/27r = 1/fc'. The transport length of a photon can 
be estimated by the absorption coefficient because 
in the considered media there is no real absorp- 
tion and attenuation of the coherent component of 
light connects only with exit of photons from the 
corresponding mode, i.e. with scattering. Taking 
into account that the amplitude of the field de- 
creases twice slower the light intensity, the ratio of 
transport length to wave length can be estimated 
as k' /2k" . Fig. 3 shows the spectral dependence 
of this ratio for clouds with different densities. 



In the spectrum of the real and imaginary parts of 





FIG. 3: Spectral dependence of the lofFe-Regel param- 
eter for atomic ensembles with different densities. 




It is seen that even for a density rt = 0.1 there is 
a region of frequencies where the considered ratio 
is less than unity. As density increases, the width 
of the corresponding spectral region also increases. 
The minimal value of the ratio fc'/2fc", however, 
decreases very slowly and shifts into the blue wing. 
It is also noticeably greater than that predicted for 
the case of independent scatterers (see for compar- 
ison Fig. 3 and Fig. 4 in [2^. Here we see directly 
the influence of the resonant dipole-dipole interac- 
tion. Increasing the density results in a decreasing 
of the portion of atoms which interact effectively 
with the light at a given frequency. 

B. Dielectric constant and atomic 
susceptibility of dense cold atomic gases 

With knowledge of the complex refractive index 
we can also calculate the dielectric constant. The 
latter can be found by the following relations 

e' = Reie) = {k'^ - k"^)/kl; 

e" = Im,{e) = 2k'k"/kl. (3.2) 



FIG. 4: Imaginary (a) and real (b) parts of the di- 
electric constant for atomic ensembles with different 
densities. 

the dielectric constant (Fig. 4) we see all regulari- 
ties which were previously observed in the spec- 
trum of the complex refractive index. But there 
is one important additional difference. For dense 
cold atomic media the real part of e can be neg- 
ative at some frequencies. The dipole dynamics 
is in phase opposition with the driving field. At 
the considered densities however, and in the cor- 
responding spectral area the imaginary part e" is 
not negligible and the electromagnetic field keeps 
wave nature. 

The complex refractive index as well as the di- 
electric constant are used for macroscopic descrip- 
tion of the light in media. One of the main char- 
acteristics in the microscopic approach is a single 
atom polarizability a which is the proportional- 
ity factor between averaged dipole moment of an 
atom and averaged strength of electric field acting 
on it. The difference between the free atom polar- 
izability and the polarizability in the medium per- 
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mits us to analyze the mutual influence of atoms 
in the medium. The key point in the calcula- 
tion of the polarizability is the idea of an effec- 
tive field acting on the atoms and its distinction 
from the mean field. In this work we will use the 
well-known Lorentz-Lorenz formula connecting the 
mean atomic polarizability and the dielectric con- 
stant \22i. 



a 



3 e- 1 
4?™ £ + 2 ' 



(3.3) 



Substituting the known dielectric constant in this 
equation give us the real and imaginary parts of the 
polarizability. The results of corresponding calcu- 
lation are shown in Fig. 5. The main result here is 
the essential decreasing of the polarizability as the 
density increases. At higher densities the collective 
atomic states are distributed over a wider region of 
frequency and polarizability at a given frequency 
per one atom is smaller. 
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FIG. 5: Imaginary (a) and real (b) parts of the atomic 
polarizability in ensembles with different densities. 



To conclude this section, let consider the quality 
of the derived dielectric constant and its applica- 
tion in a macroscopic approach. We compare two 



different results for the total cross section of light 
scattering from a homogeneous sphere. The first 
result is obtained in [l^ by means of a microscopic 
calculation. The second one is the cross section cal- 
culated in the framework of the well known Debye- 
Mie model with the permittivity from this paper. 
The result of comparison is shown in the Fig. 6. 
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FIG. 6: Spectrum of the total cross section of light 
scattering from a spherical cloud with radius R = 15. 
The atomic density is n = 0.2. The first curve is the 
results of consistent microscopic [19] calculations. The 
second curve is calculated on the basis of Debye-Mie 
model with the permittivity obtained in this paper. 



The quantitative difference between these two 
results does not exceed a few percents. It is very 
good agreement especially taking into account the 
approximate numerical determination of the per- 
mittivity. We have also noticed that Debye-Mie 
model is exactly valid for a homogeneous sphere 
whereas our atomic cloud has boundary regions 
with different local permittivity and hence is not 
completely uniform. 



IV. CONCLUSION 

In the present paper we consider the influence of 
the resonant dipole-dipole interatomic interaction 
in dense atomic clouds on their optical properties. 
Dispersion of the permittivity and atomic polariz- 
ability are determined under different conditions. 
Atomic clouds with densities up to n = 0.5 are con- 
sidered. It is observed that for a dense cloud the 
real part of the dielectric constant can be negative. 

The expression for the dielectric constant found 
here was used for calculation of the resulting spec- 
trum with that of a previous self-consistent ap- 
proach [l^. This comparison was restricted to 
the case of ensembles containing several thousand 
atoms, but good agreement allows us to use the 
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obtained permittivity for macroscopic calculations 
in cases when the microscopic approach can not 
be utilized because of technical difficulties, as in- 
dicated in the introduction to this paper. 

In this work we also determine the spectral re- 
gions, for each considered atomic density, for which 
the mean free path is smaller than the wavelength, 
i.e. we specify conditions when the loffe-Regel cri- 
terium for strong localization of light in cold dense 
atomic gases is satisfied. However more definitive 
conclusions about the possibilities of strong local- 
ization need additional study of its direct manifes- 
tation in, for instance, the distribution of fluctu- 
ations of the transmitted light intensity or in the 
afterglow delay. 

Finally, all calculations in this work were made 
under the assumption of motionless atoms but for 
non degenerate gases. In our opinion the devel- 
oped approach is applicable to the case of quasi- 
resonant compressible dipole [l^ or quasistatic 
electric dipole traps f^lj. The typical temperature 
of 30-100 mK achieved for the dense Rb cloud in 
such traps (see for example [2lj) is large enough to 
ignore all effects of degeneracy which can strongly 
affect light scattering from quantum gases [Hi . On 
the other hand, the atomic velocity is sufficiently 
small here to neglect the Doppler shift (it is sev- 
eral times smaller than the natural width of the 



excited atomic levels) and to allow us to consider 
the dipole-dipole interaction as resonant. Averag- 
ing over all possible random position of the mo- 
tionless atoms in our model allows us to take into 
account the residual motion of real atoms in the 
traps. 

It seems important to further generalize the de- 
veloped approach to the cases when atomic mo- 
tion plays a more significant role, for example to 
the case of hot gases. Such generalization is im- 
portant for a wider range problems of precision 
spectroscopy, particularly spectroscopy of selective 
reflection from the boundary of a dielectric-dense 
atomic gas [44-50]. In this case, however, the 
dipole-dipole interaction loses its resonant behav- 
ior and collisional broadening should be taken into 
account, along with essential Doppler effect. 
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